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Cell size regulation in bacteria 

Ariel Amir 1 

1 Department of Physics, Harvard University, Cambridge, MA 02138, USA 

Various bacteria such as the canonical gram negative Escherichia coli or the well-studied gram 
positive Bacillus subtilis divide symmetrically after they approximately double their volume. Their 
size at division is not constant, but is typically distributed over a narrow range. Here, we propose 
an analytically tractable model for cell size control, and calculate the cell size and inter-division 
time distributions, and the correlations between these variables. We suggest ways of extracting 
the model parameters from experimental data, and show that existing data for E. coli supports 
partial size control, and a particular explanation: a cell attempts to add a constant volume from the 
time of initiation of DNA replication to the next initiation event. This hypothesis accounts for the 
experimentally observed correlations between mother and daughter cells as well as the exponential 
dependence of size on growth rate. 

PACS numbers: 87.17.Ee, 87.17.Aa, 87.10.Mn, 87.18. Tt 



Microorganisms such as bacteria come in a diverse set 
of shapes and sizes. Nonetheless, individual strains have 
remarkably reproducible shapes, and a narrow distribu- 
tion of sizes PHI]. Many bacteria, such as E. coli, are 
rod-shaped, and during their exponential growth phase 
they elongate while maintaining a constant diameter. Af- 
ter approximately doubling their length (as well as mass 
and volume), and completing DNA replication for their 
offspring, they divide symmetrically into two approxi- 
mately identical daughter cells. In spite of decades of 
research, we still do not have a good understanding of 
how cells regulate their shape, both mechanically (i.e., 
what is the biophysical feedback necessary to achieve a 
rod-shape cell? [5]) and dimensionally: the coefficient of 
variation (standard deviatiommean, CV) can be as low 
as 0.1 for bacteria [2]. Bacteria are also remarkable in 
their ability to have a generation time that is shorter 
than the time it takes them to replicate DNA: doubling 
time Td for E. coli in rich media at 37°C is about 20 mins, 
while T r R3 60 mins are needed from initiation of DNA 
replication to cell division. This apparent paradox is ex- 
plained by the existence of multiple replication forks: in 
these situations, a cell will already start replicating DNA 
for its 4 granddaughters (or 8 great-granddaughters), in 
order for the replication to complete in time. 

Many models for cell size regulation exist in the lit- 
erature [TJ [21 IMI4], Different strategies will yield par- 
ticular cell size and inter-division time distributions, as 
well as distinct correlations. Hence, it is important to 
understand the connection between different regulation 
models and the resulting distributions and correlations. 
Moreover, there are two seemingly contradictory results 
in the literature: the first is the model by Donachie [15] . 
which shows that the measured exponential dependence 
of bacterial size on growth rate [TQ is consistent with 
initiation of DNA replication at a constant, growth-rate- 
independent volume per replication fork - suggesting a 
mechanistic picture in which a cell "knows" of its size and 
initiates replication when reaching a critical one. This 



model would imply that size at birth and division would 
not be correlated: since the time from initiation to di- 
vision is constant [17j . the size at division will be inde- 
pendent of the size at birth. However, experiments show 
that there are strong correlations between the two |18j . 

We will show here how these two results can be ele- 
gantly reconciled within a minimal model, which will be 
analytically tractable. We will suggest a mathematical 
framework which is able to capture and extend several 
existing models, and will use it to analyze the correla- 
tions and cell size distributions. We shall show that the 
aforementioned experimental data for E. coli supports 
a mechanism of cell size regulation in which the cell at- 
tempts to add a constant volume from the event of initia- 
tion of DNA replication to the next initiation event |19j . 
This model will be consistent with the results discussed 
in Ref. [TS] , predicting an exponential dependence of cell 
size on growth rate, but will also quantitatively account 
for the positive correlations between size at birth and di- 
vision [18] and negative correlations between size at birth 
and inter-division time [14] . We will show that for size- 
additive noise the size distribution is Gaussian, while for 
time-additive (i.e., size- multiplicative) noise the resulting 
size distribution is log-normal - and hence right-skewed. 
As shown in the Supplementary Information (SI), exper- 
imentally measured distributions are indeed skewed, and 
for this reason we focus on the analysis of time-additive 
noise in the main text and defer the size-additive case to 
the SI. The standard deviations of both size and inter- 
division time distributions are controlled by a single pa- 
rameter. 

The tools which we shall use will parallel those used 
when solving problems in statistical mechanics, in par- 
ticular those involving Langevin equations [50]. Mul- 
tiplicative noise and the log-normal distributions which 
emerge from our model also occur in other problems in 
physics, such as relaxations in glasses [21] and the mod- 
elling of financial markets J22j[23]. However, in contrast 
to most physical systems, negative feedback (i.e., control) 
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is a necessary feature of biological systems, including the 
problem studied here. 

Exponential growth of a single cell and regulation 
models.- The question of the mode of growth of a sin- 
gle bacterium has been a long standing problem, with 
linear and exponential growth the most common mod- 
els considered [TJ [5J [23] . Recent experiments show that 
individual cell volume grows exponentially, for various 
bacterial strains [3J l25ll2 7j . In fact, if cells grow at a 
rate that is proportional to the amount of protein they 
contain 28, 29:, as long as the protein concentration is 
constant, the cells will grow exponentially in mass and 
volume. We shall assume exponential growth of volume 
throughout this paper, v(t) cx 2*/ Td , and neglect fluctu- 
ations in the growth rate. Furthermore, we will assume 
that cells divide precisely in half since experimental re- 
sults [3D] show that division occurs at the mid-cell to an 
excellent approximation. 

Cells need a feedback mechanism that will control their 
size distribution. If cells grew for a constant time t = Td, 
random fluctuations in the timing would make the size of 
the cells at division, Vd, perform a random walk on the 
volume axis, and thus this mechanism does not control 
size. Another regulatory strategy is that of division at 
a critical mass, or of initiation of DNA replication at a 
critical size. These ideas are prevalent in the literature 
[T] [2] , but we will show that existing experimental data 
for E. coli argues against them. We shall consider the 
following class of models: upon being born at a size v b , 
the cell would attempt to grow for a time T(v b ) such 
that its final volume at division is Vd = If the 

function /(i>t>) = const, we are back to the critical mass 
model. The constant time model can also be cast in this 
language: since the growth is exponential, attempting to 
grow for a prescribed, constant time Td is the same as 
having f(v b ) = 2vb- Another important model that has 
been suggested is the so-called "incremental model", in 
which the cell attempts to add a constant volume vq to 
its newborn size [31]. In this case: 



f(vb) = v b + A. 



(1) 



In the following, we suggest a method through which 
an arbitrary regulatory model described by a function 
f(vb) can be approximately solved, i.e., we can find all 
the involved distributions and correlation coefficients an- 
alytically, finding excellent agreement with the numeri- 
cally exact solutions. We also provide methods to extract 
the model parameters from experimental data. 

The model.- We assume that the cell attempts to di- 
vide at a volume Vd — as previously explained, by 
attempting to grow for the appropriate amount of time 
t a which is a function of Vd- We assume that to this time 
is added a random noise t n , which we assume to be Gaus- 
sian. The magnitude of this noise will dictate the width 
of the resulting size and inter-division time distributions. 
Thus we have: 



^growth =t a +t n =T d log 2 [f(v b )/v b ] + t n , 



(2) 



with t n assumed to be a random variable with: P(t n ) = 

, 1 = e 2 "t . The model is similar to that discussed in 

Ref. [32] , where the molecular mechanisms leading to the 
noise in budding yeast are studied. 

We will calculate the inter-division time and volume 
distributions. The key insight is that for noise that is 
not too large (equivalent to size distributions which are 
not too broad, i.e., with a small CV), it is the behavior 
of f(i>b) around the average newborn size v that is the 
most important. Therefore we can Taylor expand f(vb) 
around v: 



f(v b ) «/(«) + f'(v)(v b -v). 



(3) 



As an example, the incremental model has /' («) = 1 and 
V = A, while the critical size model has f'(v) = 0. 

Any two models that agree to lowest order, will result 
in similar distributions - provided the noise is not too 
large. We therefore choose to solve an equivalent model, 
that will be amenable to analytic treatment, and that 
can be viewed as an interpolation between the critical 
size model and the constant doubling time model. We 
choose: 



t a = Td[l + a]og 2 (v 0 /v b )], 



(4) 



which corresponds to the regulatory function: f(v b ) = 



v 0 2 t ° = 2vl 



The case a — 0 corresponds to con- 



stant doubling time model (/'(t>o) — 2), while a = 1 
corresponds to the critical size model (/'(^o) = 0). Im- 
portantly, for a — 1/2 we have /'(i>o) = 1, as does the 
incremental model: hence, using a target function like 
this gives results close to a perfect realization of the in- 
cremental mode. We shall show that the paramter vq in 
Eq. (14 1 will be very close to the average newborn cell 



size v. 

Solution of size and inter- division time distributions. - 
We shall consider the case of symmetric division, relevant 
for many bacteria. For a newborn size v b , we have for the 
next newborn volume: vl 
Therefore: 



..new 

J b 



v 0 v b z 



log 2 «*>o) = (1 - a) \og 2 (v b /vo) + t n /r d , (5) 



From stationarity of the stochastic process we know 
that P{v^ ew ) = P(vb). Since t n is a Gaussian variable, 
we find that log 2 (u&) is also a Gaussian variable, and 
hence P{v b ) would be a log-normal distribution. If we 
denote the variance of log 2 (?;f,/i>o) by a^, we have cr^ = 
cr^fl — a) 2 + therefore the newborn size distribution 
is: 



Downloaded from http://biorxiv.org/on September 18, 2014 



3 



[l°82(V"0)r 
2^ 



P(v b ) 



2ir ln(2)cr t 



Vb 



with 



rj a {2 - a) ■ 



(6) 



(7) 



Note that the average cell size is v — voe ln ( 2 )°"/ 2 ; for 
realistic values of o~ v it will only be a few percent larger 
than Vq. Similarly, the standard deviation of the size 
distribution will be approximately a s ~ ln(2)<7„Vo, and 
the coefficient of variation is thus: vcv ~ ln(2)<T t ,. The 
skewness of the distribution is positive: 71 3\n(2)(r v , 
and provides a useful test of the assumption of a time- 
additive rather than size-additive noise, as we elaborate 
on in the SI. 

We can now find the distribution of division times us- 
ing: trf = t a + t n . Since Vf, depends only on the noise of 
previous generations, t a is independent of t n , and since 
log 2 (vt)/t>o) and t n are Gaussian variables, the resulting 
inter-division time distribution is also Gaussian, and has 
a variance given by: 



Var[t a 



222 
r^a <7„, 



2- a 



(8) 



In the case a — > 0, we find that a v diverges (an ex- 
tremely broad distribution of newborn sizes), but the 
inter-division time distribution is narrow: Var[i<i] ~~ ► 
ax 2 , as should clearly be the case since there is no size 
feedback mechanism in this case. Note that for any pos- 
itive value of a > 0 there will be a stationary size distri- 
bution, but for a = 0 there is no stationary distribution. 



From Eq. (|16|) we find that the CV of the distribution 

'—. It is 



of inter-division times is given by: tcv 
instructive to consider the dimensionless quantity: 



Trf 



7 = vcv /tcv 



ln(2) 
V2a' 



(9) 



By constructing 7 from the experimental distributions 
we can extract the value of a and find the form of size 
regulation utilized by the organism, if the division is sym- 
metric. Later we shall show an additional, independent 
way of extracting a, which will be more robust against 
measurement noise since it will rely on correlations rather 
than the distribution widths. 

Fig. [l] compares the numerically obtained size distri- 
bution for various values of a and the incremental model, 
with the result of Eq. ([6]), finding excellent agreement. In 
the SI we extend this comparison to various noise magni- 
tudes. Our model captures the numerically exact solution 



£b o Numerics, a = 0.1 

Numerics, a = 0.5 
■ Numerics, a = 1 
1 Incivnu'iita] imxU'l 




Vb (a.u.) 




t (mins) 



FIG. 1. Comparison between the analytical results of the 
model for varying values of a (Eqs. ( 6)16 1), and numerics. 
Choosing a — 1/2 provides an excellent approximation for 
the incremental model, as the effective size regulation of the 
two models agrees to lowest order. The parameters of the 
model are chosen according to their realistic values for E. coli 
growing at 37°: doubling time is Td = 20 mins and ar/rd = 
0.2 [2]. For each case, the numerical distribution is extracted 
from a sequence of 10 7 divisions. 



very well and Eq. ( 14 ) provides a useful tool to capture 



a generic division strategy characterized by an arbitrary 
function /(«&). 

Extracting the parameters from experiments.- Within 
the class of models proposed here, the value of a can be 
obtained by considering the correlations between size at 
birth and size at division. In the SI we show that the 
Pearson correlation coefficient between size at birth and 
division (equivalent to the mother: daughter correlation 
coefficient) equals: 



Csal-a. 



(10) 



In particular, for the incremental model the correlation 
coefficient between mother and daughter cells should be 
0.5. 

Upon fixing the value of a, a single parameter, ctt, 
will determine the distributions of both size and division 
time, and the calculations performed here would allow 
one to scale both distributions using this single parame- 
ter. For the time-additive noise analyzed here, the model 
predicts an approximately log-normal newborn size dis- 
tribution, given by Eq. ^ , and a Gaussian inter-division 
time distribution, given by Eq. (16), whose standard de- 



viation is larger than <tt- In the SI we show that for size- 
additive noise, one obtains a skewed time-distribution 
but a Gaussian size distribution - in contrast to what 
is observed experimentally |13) . Therefore, observing the 
distribution shape provides useful information regarding 
the source of the noise. Further experiments are needed 
to elucidate the molecular source of this multiplicative 
noise. 

Cell size control in E. coli.- Experimentally, various 
correlation coefficients were measured for E. coli at slow 
growth conditions in Ref. [18] . using the membrane elu- 
tion technique. The correlation coefficient between new- 
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born cell and daughter cell size was found to be C = 0.55, 
close to the theoretical 1/2 value expected for the in- 
cremental model. There was a strong correlation (0.8) 
between size at initiation of DNA replication and size 
at division, as we would expect from the assumption 
of exponential growth and that the time from initiation 
to division is constant [T7]. Yet these observations ap- 
pear to be in direct contradiction to the idea that ini- 
tiation occurs at a critical size [15]. The key point is 
that Donachie's analysis [T5] shows that there is a criti- 
cal size for initiation of DNA replication (independent of 
growth rate) , on average. It is only from the fluctuations 
(i.e., correlations) that one can understand whether the 
underlying regulatory mechanism utilizes a critical size 
or integrates volume - as we shall propose is the case. 
Ref. [IH] gives a simple biophysical implementation of the 
incremental model, which will reconcile these seemingly 
contradictory results and will realize a particular case of 
the class of model we proposed here: in this model, a 
protein A is forced to have a growth-rate-independent 
density throughout the cell using a negative feedback in 
its regulation, and a second protein B is produced when- 
ever A is. In this way when cell volume grows (and only 
then), more A and B proteins are generated in an amount 
proportional to the change in volume. The hypothesis is 
that B proteins localize at their potential initiation site 
(which we will assume to be one of the replication ori- 
gins), and only when their total number at each origin 
reaches a critical value does initiation of DNA replication 
occur, after which B is degraded. Note that two types of 
proteins are necessary, since in order to measure volume 
differences A must be spread throughout the cell, while 
B has to localize to measure an absolute number (rather 
than concentration). See Ref. [T!5] for further details. 

In the SI we show that this model reduces to the in- 
cremental model, albeit with an effective A (see Eq. 0) 
which strongly depends on the growth rate: 

A = A2 T "/ Td . (II) 

According to our results the average cell size will be 
A - in agreement with the experimental results seeing 
precisely this exponential dependence of bacterial size on 
growth rate, with T r the exponent |16] . This model nat- 
urally accounts for the "quantization" of the cell critical 
mass at initiation at different growth rates |15j . without 
necessitating the measurement of an absolute mass or 
volume. Moreover, it is plausible that the source of noise 
in adding the incremental volume will be due to "molecu- 
lar noise" (number fluctuations of protein B), and would 
therefore be weakly dependent on growth rate. The same 
calculation which leads to Eq. |2l] ) would suggest that 
o~t (the noise standard deviation) should depend on the 
growth rate in the same exponential way as A. This im- 
plies that the CV of size distributions should be weakly 
dependent on growth rate (see Eq. (|7|), an expectation 
supported by Ref. [55] , 



Thus, we have shown that using our calculations and 
the interpretation in terms of the incremental model ex- 
plains various experimental results. In fact, the model 
also makes precise predictions with regards to additional 
correlations: for example, it is possible to show that for 
the incremental model the size correlation coefficient be- 
tween cells N generations apart is 2~ N . Similarly, the 
model predicts a negative correlation of -1/2 between the 
size at birth and the inter-division time, see the SI for 
further details. This correlation coefficient was recently 
analyzed by Robert et al. [H], using data from two differ- 
ent experimental systems (with different growth rates), 
finding a correlation coefficient of -0.5 in both cases: ex- 
actly as predicted by our model. This gives a particularly 
simple and transparent interpretation to their analysis, 
and provide additional, strong support for the incremen- 
tal model. Refs. [T3J and [33] find similar negative cor- 
relations between newborn size and inter-division size, 
supporting our conclusion. 

All of these provide additional support for the rele- 
vance of this model to cell size control in E. coli, and 
most likely to other organisms as well. It is possible, how- 
ever, that alternative biophysical mechanisms may lead 
to the same correlations and size dependencies calculated 
here, and for this reason finding the underlying biological 
mechanism is important; In recent years, dnaA has been 
shown to have properties reminiscent of the biophysical 
model described here [5], where its active and inactive 
forms correspond to the roles of proteins A and B above 
- see the SI for further details. 

Discussion.- In this work we suggested a phenomeno- 
logical model which is able to describe partial size control 
within a broad class of control strategies, and interpolate 
between the case of constant time to division and division 
at a critical size, for both size-additive and time-additive 
noise. We are able to analytically calculate the size and 
inter-division time distributions for the case of symmet- 
ric division, relevant to various bacteria. For E. coli, we 
have shown that a simple biophysical model in which a 
constant volume is added from consequent events of ini- 
tiation of DNA replication predicts: I) Cell size depends 
exponentially on growth rate. 2) Cell size distributions 
are approximately log-normal. 3) The correlation coeffi- 
cient between size at birth and division is approximately 
1/2. 4) The correlation coefficient between size at birth 
and time to division is approximately -1/2. 5) The ra- 
tio of the CV of size and inter-division time distributions 
is approximately ln(2). The simplicity of a biophysical 
model which implements this idea [19] suggests that this 
may be a robust way of regulating cell size and coupling 
DNA replication and growth. 

This interpretation in terms of the incremental model 
suggests an outstanding puzzle: can we underpin the pre- 
cise molecular mechanism responsible for volume integra- 
tion? Can the source of the noise in inter-division times 
be elucidated? Testing this model further in other mi- 
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croorganisms may yield important insights into cell size 
regulation, and in particular, it is intriguing to see if the 
same ideas are applicable to cell size control in higher 
organisms. Recently, size distributions in other microor- 
ganisms were shown to obey simple scaling laws [35] , sug- 
gesting this to be a promising direction, and that the 
model discussed here may have a broader range of appli- 
cability. 
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CELL SIZE REGULATION IN 
MICROORGANISMS - SUPPLEMENTARY 
INFORMATION 



Calculating the correlation coefficients 

In the following we calculate the Pearson correlation 
coefficients between various variables of interest, for the 
model with an arbitrary value of a characterizing the size 
control strategy. 

Mother size '.daughter size correlation coefficient 
For a narrow size distribution (such as the biologically 
relevant one) , performing a linear regression analysis be- 
tween the size at birth and the size at division would 
not be very different than doing the regression between 
x = \og 2 {v b /vo) and y = \og 2 {vd/vo) (for the data corre- 
sponding to Fig. 1, for example, the difference between 
the two methods is less than one percent). Within our 
model, the value of the dimcnsionlcss slope j3 can be read- 
ily found from Eq. (3) of the main text: since the noise 
is uncorrelated with the random variable x, we have: 



Thus: 



C ts » E(t d \og 2 (v 0 /v b ))/cr t a v = - 



a>T d Var[\og 2 (v 0 /v b )} 

yJVar[t d ]o- v 

(15) 

with Var[td] the variance of the generation time which 
was calculated in the manuscript and found to be (Eq. 
(7) of main text): 



Var[t d ] = a T ' 



2-a' 



(16) 



where ctt is the noise magnitude. 

Similarly, the variance of the logarithm of newborn 
size, er^, was found to be (Eq. (6) of main text): 



Var[\og 2 [v Q /v b ] = a T 2 jr\ 



Combining these we find: 



a(2 - a) 



(17) 



/3 = />o) = 2(l-a). (12) 

(i.e., it is 1 for the incremental model). For symmet- 
ric divisions, the slope of the linear regression between 
a newborn cell and the size of the daughter cell immedi- 
ately after division would thus equal 1—a. The coefficient 
of correlation between the newborn cell and the newborn 
daughter cell is: 

C'md = E(xy)/a x a y = (1 - a)E(x 2 )/a x a y . (13) 

From stationarity, we know that the distribution of v b 
and of v^ ew are identical, hence a y — <r x , and: C = 1 — 
a. Therefore, for the incremental model the correlation 
coefficient between mother and daughter cells should be 
0.5. 

Newborn Size: Generation time correlation coefficient 
From Eq.(4) of the main text: 

t a = T d [l + alog 2 (v 0 /v b )], (14) 

Hence, we expect a negative correlation between 
log(newborn size) and inter-division time, with slope —a. 
Such a negative correlation was recently observed pQ. 
Let us calculate the Pearson correlation coefficient of the 
inter-division time with newborn size, not to be confused 
with the aforementioned correlation. Since the distribu- 
tion is quite narrow, this is nearly the same as the cor- 
relation coefficient between logarithm of size and inter- 
division time, following the same logic which we used for 
the evaluation of the mother size:daughter size correla- 
tion coefficient. 



Cts = -^1 ■ (18) 

In particular, it is -1/2 for the incremental model. 

Reduction of biophysical model to the incre- 
mental model 

We shall now show that the biophysical model of Ref. 
[2] , where a constant volume is added between two DNA 
replication initiation events, perfectly realizes the incre- 
mental model (corresponding to Eq. (1) of the main), 
yet with A which depends on size in a particular way: 
Let us assume that that there are 2™ replication forks 
at work (and hence 2 n+1 replication origins), and that 
initiation of DNA replication occurred at a volume for 
one of them. Protein B will be accumulated at each ori- 
gin until a critical amount is reached. This implies that 
the next initiation will occur (on average) at a volume 
v next _ _|_ 2 n+1 A)/2, where according to the above 
biophysical mechanism A is independent of the growth 
rate. According to our assumption, there will be n di- 
vision events from initiation to the end of replication. 
Since growth is exponential and we are assuming per- 
fectly symmetric divisions, if the cell volume at initia- 
tion is Vi its volume at the end of DNA replication is 
v d = Vi2 Tr / Td ~ n , regardless of when the intermediate di- 
vision events happened. Its daughter cell will be born 
with a volume v nb — v d /2, and its size at division will 
equal, by the same reasoning: 

ynext _ ^next^Tr/rd—n _ f} i 2 T '-/ T d~ n _|_ 2™ A2 T '"/ r<i_n 

(19) 
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relevant narrow distributions will be essentially equiva- 
lent, a — 0 has no size control, a — 1 realizes division 
at a constant size, and most importantly, a — 1/2 is an 
exact realization of the incremental model. 

Now we assume that the noise is added to size rather 
than time, hence: 



(23) 



with £ a Gaussian variable with standard deviation erg. 



Using Eq. (27 1 we have: 



FIG. 2. Numerical simulation compared with Eq. (6) of 
the main text, for the incremental model (corresponding to 
a — 1/2), for varying noise. The noise is added to the inter- 
division time and has a standard deviation gt- For each case, 
the numerical distribution is extracted from a sequence of 10 7 
divisions. 



Thus we have (up to the noise): 

next 1 a 

Vd = v nb + A, 

with: 

A = A2 Tr/Td . 



(20) 



(21) 



vT w = av 0 + (1 - a)v b + t (24) 
which can be rewritten as: 

(v2 ew -v 0 ) = (l-a)(v b -v 0 )+Z, (25) 

Following the same reasoning used in the main text, 
v b — vq will be a Gaussian variable with vanishing mean, 
and: 



Var[v b ] = (1 - a) 2 Var[v b ] + a 2 s , 



(26) 



Hence: 



Comparison of theory and numerics for time- 
additive noise 

In this section we elaborate on the generality of the 
results derived in the main text by showing the excellent 
agreement between the analytical form of Eqs. (6) and 
(7) of the main text and the numerics. The distributions 
were evaluated numerically by following the lineage of 
10 7 divisions for each case. 

Fig. [2] compares the analytical results for a — 1/2 with 
numerical results on the incremental model, for varying 
noise. We tested noise magnitudes ranging from the bio- 
logically relevant value of ctt/t = 0.2 to relatively large 
noise with o~t/t — 0.6. We find excellent agreement even 
for large noise - despite the fact that the incremental 
model is equivalent to the a = 1/2 analytically tractable 
model only to lowest order in (v b — v 0 ). 

Solution of model for size-additive noise 

Using the tools presented in the main text, it is 
straightforward to modify the approach for size-additive 
noise. In this case, it is convenient to use: 

f{v)=2v 0 +2{l-a){v-v 0 ). (22) 

This agrees to lowest order in (v — vq) with the defini- 
tion of a of the main text, and hence for the biologically 



Fig. [3] verifies this result, comparing the analytic re- 
sults for various values of a with numerics (including 
a = 1/2 which corresponds to the incremental model), 
all with size-additive noise with as/ v o = 0.27 (chosen to 
capture the width of the experimentally observed distri- 
bution). 

Distinguishing between size-additive and time- 
additive noise 

For the realistic biological parameters, the difference 
between the time-additive and size-additive noise is not 
dramatic, as is shown in Fig. [4] 

We also tested the size distributions resulting from the 
regulation strategy corresponding to a generic value of 
a (Eq. (3) in the main text), with a noise that has 
both a time-additive component (with standard devia- 
tion o~t) and a size- additive component (with standard 
deviation 0-5). Based on the results described above that 
an approximately Gaussian distribution occurs for a size- 
additive noise, we expect that only the time-additive 
component of the noise will significantly contribute to 
the skewness of the distribution, since a Gaussian distri- 
bution has vanishing skewness. Therefore, we calculated 
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FIG. 3. Comparison between the analytical results of a Gaus- 
sian distribution with variance given by Eq. ( 27 1 and numer- 



ics with size- additive noise, for varying values of a. Choos- 
ing a = 1/2 provides an exact realization of the incremental 
model. 



1.8 
1.6 
1.4 
1.2 
1 



Numerics, size-additive noise 
■ Gaussian, Eq. (6) 
Numerics, time-additive noise 
Log-normal, Eq. (5) of main text 




Vb (a.u. 



FIG. 4. The distribution of newborn cell volume was found 
numerically, using the incremental model for both the cases 
of size-additive and time-additive noise. In the first case 
os/v 0 = 0.21n(2) « 0.14, while in the latter a T /r d = 0.2. 
As is shown, the two cases are well-approximated by Gaus- 
sian and log-normal distribu tion s, respectively, and by the 
theory corresponding to Eq. (271 of the SI and of Eq. (6) of 
the main text. 



numerically the skewness of the newborn size distribu- 
tion: 71 = E(vb — /i) 3 /<7 3 , where /i w vq is the average 
newborn size and a is the standard deviation of the size 
distribution. 

As expected, we found that in order to have significant 
skewness, there has to be a time-additive noise compo- 
nent present, as is illustrated in Fig. [5] 

Interestingly, in the analysis of financial markets a sim- 
ilar strategy is used to test the multiplicative nature of 
the stochastic processes underlying the market fluctua- 
tions pHi]. 




n 



I 



FIG. 5. The skewness of the newborn size distribution is 
shown in the 2d map for the mixed case where both time- 
additive and size-additive noise are present, with magnitudes 
ot and as- As is shown, distributions with non- negligible 
skewness only occur in the presence of time-additive noise. 



Experimental support for time-additive noise 

The purpose of this section is to show that existing ex- 
perimental data for E. coli supports time-additive (i.e., 
multiplicative) noise rather than size-additive noise, us- 
ing the tool developed in the previous section: namely, 
the skewness of the size distribution. 

Ref. [5] measures the newborn size distribution of E. 
coli using the membrane elution technique, and shows 
that it has non-negligible positive skewness and it agrees 
very well with a log-normal distribution. 

Recently, Ref. [1 analyzed cell size fluctuations 
and distributions. Their size distributions are strongly 
skewed to the right (71 > 0.5), suggesting that the noise 
is multiplicative rather than additive to the size. A sim- 
ilar analysis was performed in Ref. [5], which further 
corroborates our conclusion. 

Although these various experiments suggest that noise 
in the size-control of bacteria is multiplicative rather than 
additive to size, a significant amount of additional exper- 
imental work is needed in order to elucidate the origin of 
the noise at the molecular level. Ref. [7] used the pow- 
erful tools of molecular biology to gain insights into pre- 
cisely such a problem in budding yeast. In this pioneering 
work, they also model the noise as time-additive rather 
than size-additive, and are able to deduce the molecu- 
lar mechanism for it and modify it in order to test their 
stochastic model. It would be highly rewarding to do the 
same for bacteria. 

DnaA as a volume integrator 

Here, we review some recent progress in our under- 
standing of the function of dnaA in initiating DNA repli- 
cation in bacteria, and show that it shares many proper- 
ties with the hypothetical model described in the main 
text. DnaA has two forms, an active, ATP-bound form, 
and an inactive ADP-bound form. It is believed that 20 
dnaA proteins in their active form are needed in order to 
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initiate DNA replication [5] [3] ■ This highly cooperative 
mechanism is parallel to the "critical number" of protein 
B in the main text (P2 in Ref. [2 ). After initiation, 
the number of active copies of dnaA quickly drops [9], 
which is equivalent to the "degradation" of protein B. 
Importantly, dnaA is known to autoregulate [10] . which 
in certain cases leads to a concentration approximately 
independent of the growth rate [11] - this was the re- 
quirement for protein A of the main text (Pi in Ref. 
[5]). As expected, inducing high levels of dnaA in E. coli 
leads to early initiation |12j . while dnaA mutants initi- 
ate at a later time [13] [TJ]. Therefore, it is possible that 
the two forms of dnaA serve together to implement the 
volume integration which is hypothetically discussed in 
Ref. [2J. Certainly, many other proteins are involved in 
the process, such as SeqA [13| . and further experiments 
are needed to better understand the biochemical and bio- 
physical mechanisms. Nevertheless, the results presented 
in the main text lead to severe constraints on possible 
models. 
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